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Abstract 

The classification of lattice equations that are integrable in the sense of higher- 
dimensional consistency is extended by allowing directed edges. We find two cases 
that are not transformable via the 'admissible transformations' to the lattice equa- 
tions in the existing classification. 

PACS: 02.30.Ik, 04.60.Nc, 05.45.Yv 

1 Introduction 

Integrable systems come in many forms. It can be argued that of these different forms 
(which include partial and ordinary differential equations, differential-difference equa- 
tions, and difference equations) lattice (or 'partial difference') equations are the most 
fundamental: they include not just one continuous integrable equation in the continuum 
limit, but entire hierarchies (see Q, |[T3]| ). For systems defined on the vertices of a two- 
dimensional quadrilateral lattice or, more generally, on quad-graphs, the definition of 
integrability is so natural and transparent that it adds further weight to the argument for 
the primacy of such systems in the pantheon of integrability. For these systems integrabil- 
ity is synonymous with higher-dimensional consistency. Integrable lattice equations are, 
in a sense, their own Lax pair (or 'zero-curvature condition') - the Lax pair follows in an 
algorithmic manner from the lattice equation (and the lattice equation follows from the 
Lax pair). The clarification of three-dimensional consistency as the definition of integra- 
bility for such systems and the exploitation of the definition for the derivation of the Lax 
pair is a recent development, lITOl . Ol . ISl (we refer to these papers for further details). 

The importance of such systems established, it becomes desirable to have a classifi- 
cation, and a notion of which transformations can be performed on a lattice equation such 
that it still belongs to the same position in the classification. These issues were coirfronted 
for the first time in the landmark paper m. 
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In |[T1, Adler, Bobenko and Suris ("ABS") classified 'integrable' equations on quadri- 
lateral lattices under certain symmetry assumptions (where integrability is equated with 
thi^ee-dimensional consistency) and specified the admissible transformations that may be 
performed on these equations such that they are still considered to be the same system 
(that is, take the same place in the classification). 

The ABS classification has already received further study: reductions of the 'Q' list to 
second-order mappings was performed by Joshi et al. in ||6l, in Hietarinta studied the 
necessity of the 'tetrahedron property' assumption, and symmetries of the lattice equa- 
tions in the classification were investigated by Rasin and Hydon in \\\2l . The classification 
has influenced other papers to varying degrees. 

In the present paper a symmetry assumption of ABS is generalized in a natural way. 
This extends the classification by the addition of the two integrable lattice equations 

Q={x-y){v-u) + {a + ^){u-v) + {^- a){x -y)+ X{a) - ?i(p), 

Q = -n{-^)xu - ^l{^)vy + ^i{-a)xv + ^i{a)uy, 

where X is any even function, and n is any function that is not even. We call these 
equations Bl and B2 respectively. They ai^e inequivalent to the existing members of 
the classification (they cannot be transformed to existing members of the classification 
by 'admissible transformations'). The result here is that a generalization of one of the 
symmetry assumptions results in a new classification, consisting of the lattice equations 
in the original ABS classification, plus the addition of two lattice equations that obey the 
new symmetry (but not the more restrictive form stipulated by ABS). 

In section |2] the assumptions made by ABS in HI are described in detail. In section 
[3] the change to the symmetry assumptions is given; the modification to the derivation of 
ABS due to this change is given in section HI In section [5] the new lattice equations are 
derived, and they are discussed in section [6l 

2 A review of the ABS assumptions and derivation 

The equations on quadrilaterals considered by ABS are of the type 

Q{x,Xi,xj,Xij;ai,aj) =0, (1) 

where the field variables x,Xi,Xj,Xij € C are attached to the vertices of the quadrilateral 
(subscripts / and j denoting steps along the lattice in directions / and j respectively) and 
'lattice parameters' a, and are associated to edges running in the / and j directions 
respectively (see fig. 1). 
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Lattice equations of the form ([T]) are classed as integrable if they have the property 
of three-dimensional consistency. To understand what 'three-dimensional consistency' 
means, consider the lattice equation ([Hi to be defined on the faces of a three-dimensional 
combinatorial cube, where parallel edges cany the same lattice parameter. For instance, 
labelling the three orthogonal directions (1,2,3), edges in the 1, 2, and 3 -direction have 
the lattice parameter ai, aa, and a3 attached, respectively. Considering the initial data 
x,xi,X2,xt, to be given, three-dimensional consistency states that the three different ways 
of deducing xuj all give the same result (see fig. 2). 



Fig 2. Three-dimensional consistency. Given initial values at the black dots, x,xi ,X2,X3, 
the three possible ways to calculate x\23 all give the same result. 

For notational reasons, when considering only one plane (say the (1,2) plane) we will 
denote tti = a, tti = p, xi =u,X2= v, x\2 = y. 

The assumptions under which ABS solved the classification problem of such equa- 
tions are listed below. 
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1] Linearity. 

Q{x,u,v,y;a,^) is affine linear in each argument: 

Q{x, u,v,y; a, P) = aixuvy + . . . + aie 

where {«,} depend on a, p. 
2] Symmetry. 

Q{x,u,v,y;a, P) is invariant under the group of square symmetries: 

Q{x,u,v,y;a,^) = eQ{x,v,u,y;^,a) = aQ{u,x,y,v;a,^) (2) 

s,a = ±1. 

3] Tetrahedron property. 

The function 

■^123 = zix,xi,x2,xy,ai,a2,a3), 

obtained from three-dimensional consistency, does not depend on the variable x, that is, 
z. = 0. 

The following transformations are assumed to identify equivalence classes: 

• Action on all field variables {x,} by one and the same (independent of lattice pa- 
rameter) Mobius transformation. 

• Simultaneous point change of all parameters a i-^ <l)(cx). 

These transformations do not violate the three assumptions listed above. 

That these are the only admissible transformations is a key point of the classification. 
It is stated in ^ that the 'A' hst is obtainable from the 'Q' list, and the lower 'Q' and 'H' 
lattice equations (Ql to Q3, and HI to H2, respectively) are obtainable from the top equa- 
tion in their lists (Q4 and H3, respectively) by wider classes of transformations. However, 
it is stressed that these transformations are outside of the admissible transformations that 
the classification rests on. (The issue of transformations of integrable lattice equations 
was discussed again in Q, where, again, transformations on sublattices are observed to 
sometimes preserve and sometimes conflict with higher-dimensional consistency.) 

Under the ABS symmetry assumptions there are two possible forms of the lattice 
equation, corresponding to the two possible values of a in equation (O. 
a = l: 

Q = aixuvy + a\ [xuv + uvy + vyx + yxu) + (xy + uv) 
+a2{xu + vy) + a2{xv + uy) + a-i{x + u + v + y) + 04. (3) 

a = -1: 

Q = ai {xuv + uvy — vyx — yxu) + a2 {xy — uv) + aj{x — u — v + y). (4) 

Under the ABS assumptions the a = — 1 case was shown to be empty. By removing 
the assumption of the tetrahedron property, Hietarinta was able to derive an integrable 
lattice equation with this a = — 1 form; however, it was later shown to be linearizable by 
Ramani, Joshi et al. lITTI . 
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3 A change of assumption 

In this paper, assumption 2 (the symmetry assumption) is modified to allow directed 
edges. (A notion of directed edges has appeared at least once before in the literature, 
HI, although the labelling of the quadrilateral was slightly different to that given here.) 
We change assumption 2 to 

2' ] Symmetry. 
Q{x,u,v,y;a,^) is invariant under the symmetries: 

Q{x,u,v,y;a,^) =eQ{x,v,u,y;^,a) = aQ{u,x,y,v;-a,^) (5) 

s,a = ±1. 

The change here extends the reflection symmetry of ABS, {x,u,v,y) i— > {u,x,y,v), in 
a natural way to include 'directed edges', in the sense that under this reflection transfor- 
mation we also have a i-^ — a. This is a generalization of (as opposed to an alternative 
to) the ABS classification. The ABS classification consists of the set ofintegrable lattice 
equations with these symmetries where the lattice parameters {a,} appear in evenfiinc- 
tions. For example, take the members of the ABS classification as presented in 111, make 
the replacement a,- 1-^ aj, for all /, and they have this new symmetry. 

Under the new symmetry assumption ^ we can write 

Q=A xuvy + Bxuv + C uvy + D vyx + E yxu + F xy + Guv + Hxu + Ivy + J xv + Kuy 

+ Lx + Mu + Nv + Oy + P, (6) 

where A = A(a,P), etc., and the symmetry conditions impose the following conditions 
on the coefficients: 

A(-a,p) = aA(a,p) , B(-a,p) = a£'(a,|3) , C(-a,p) = aD(a,p) 

F(-a,p) = aG(a,p) , //(-a,p) = a//(a,p) , /(-a,p) = a/(a,p) 

7(_a,p) = a7i:(a,p) , L(-a,p) = aM(a,p) , N{-a,^) = aO(a,p) 
P(-a,p) = aP(a,p) 

and 

A((3,a) = £A(a,p) , S(p,a) = £B(a,p) , C(p,a) = £C(a,p) 

D(p,a) = ££■(«, p) , f(p,a) = eF(a,p) , G(p,a) = £G(a,p) 

//(p,a) = £7(a,p) , /(p,a) = £/«:(a,p) , L(p,a) = £L(a,p) 

M(p,a) = £Ar(a,p) , 0(p,a) = £0(a,p) , P(p,a) = £P(a,p). 

To preserve this symmetry the natural transformations under which equations are 
identified are now: 

• Action on all field variables {xi} by one and the same (independent of lattice pa- 
rameter) Mobius transformation. 

• Simultaneous point change of all parameters a i— > (])(a) that preserves the parity of 
the lattice equation coefficients as functions of the parameters. 
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These are our admissible transformations. When the lattice coefficients depend on the 
lattice parameters only as even functions (i.e., the lattice equation has 'undirected edges') 
then the class of simultaneous point changes of all parameters is the same as the ABS 
case; however, if the lattice coefficients depend on the lattice parameters in ways other 
than as even functions, the transformation of parameters must obey (|)(— oc) = —({)(«), or 
the reflection symmetry will be lost. Any simultaneous point change of all parameters 
will preserve integrability but may break the new symmetry (so that the lattice equation 
may have neither symmetry 2 or 2'). Note again that if the lattice equation coefficients are 
even functions of a (undirected edges) then the form of the lattice equation, Q, becomes 
that of ABS ©, ©. 

The fractional multi-linear expressions that arise when solving the lattice equation (in 
this case (O) for a particular field variable are an important ingredient of the algorithmic 
way of obtaining the Lax matrices from the equation. Lattice equations that do not lead 
to rational expressions (for instance equations without multi-linear terms) are discounted 
from the definition of integrability on account of being trivial. Lattice equations that 
do not contain lattice parameters are also excluded from the definition of integrability 
(note that in this case the Lax matrices that follow from the algorithmic construction, 
guaranteed by three-dimensional consistency, will not have a spectral pai^ameter). 



4 Derivation: Analysis 

We now follow the first part of the classification derivation in [H (the 'analysis') making 
changes where necessary for the different symmetry condition. 

The 'analysis' of the classification of ABS proceeds by using the assumptions to 
deduce a parameter-less discriminant associated with the quad-graph equation. Proofs 
identical to those of lUl are omitted. 

In m it is shown that 

g{x,xi;ai,a2) g{x,x2;a2,a3) g{x,x3;a3,ai) 
= -g{x,xi;aua3)g{x,x2;a2,ai)g{x,x3;a3,a2), (7) 

where g(x, u; a, P) is a biquadratic polynomial in x and u, defined by either of the formulas 

g{x,u;a,^) = QQyv-QyQv (8) 

g{x,vM = QQyu-QyQu (9) 

where Q = Q{x,u,v,y;a,^). We now have a slight change from ABS, 'the polynomial g 
is symmetric g{x,u;a,^) = g{u,x;a,^y (HI, page 521, equation (19)) becomes, due to 
the different symmetry, 

gix,u;a,^) =g{u,x;-a,^). (10) 

As in m, one may prove the following lemma. 

Lemma. The discriminants of the polynomials g = g{x,u;a, P) and g = g{x,v;^,a), 
considered as quadratic polynomials in u and v, respectively, coincide: 

gl - ^gguU =gl- 2gg„y . (11) 
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Equation (|7]l implies properties of the polynomial g. The following proposition is the 
same as proposition 5 of HI (the argument of the proof is similar; however, the differences 
give us additional information on k{a, P)). 

Proposition. The polynomial g{x,u;a, p) may be represented as 

g{x,u;a,^) =k{a,^)h{x,u;a), (12) 

where the factor k is antisymmetric, 

^:(P,a) = -^(a,p), (13) 

and the coefficients of the polynomial h{x, u; a) depend on the parameter a in such a way 
that its discriminant 

r{x)=hl-2hKu (14) 

does not depend on a. 

Proof. Equation ([7]) implies, by a sepai^ation of variables argument, that 

g(x,xi;ai,a2) /(x;ai,a2) 
g{x,xi;ai,a3) ' f{x;ai,aj) 

and, hence, by the new symmetry 2' and a further separation of variables argument 

g{x,xi;ai,a2) k{ai,a2) 



g(x,xi;ai,a3) k{ai,a^) 

as in m, but with the additional restriction that 

k{a\,a2) _ k{-ai,a2) 
k{ai,ai,) k{-ai,a3)' 



(15) 



(16) 



Then, continuing as in HI, equation ^ with (fTSl ) implies the antisymmetry of k, i.e., 
k{a,^) = -k{^,a), and (fTSl ) implies (fT2l ). Equation (fTTI ) implies 

h^^ — Ihhuu = h^, — 2hhvy , h = h{x,u;(x) , h = h{x,v;^) 

and therefore r(x) :=h^ — Ihh^u does not depend on a. □ 
By a separation of variables argument, equation (fT6] ) gives 

^(ai,a2) _ x(ai) ^^^-^ 



k{-aua2) x(-ai)' 

However, from (fTSl ). ^(a, P) is only defined up to a factor that is a function of a. Hence 
we may take ^(a, P) to be defined such that 

^(ai,a2) = ^(-ai,a2)- (18) 
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The following lemma is new, and is a key fact leading to the derivation of the new 
lattice equations. 

Lemma. The biquadratic h is symmetric in x and u, 



h{x,u,a) = b(jX u +bi{x u + xu ) + b^{x +u ) + b5Xu + b(,{x + u) + bs, (19) 

and {bi} are all even functions of a. 

Proof. The factorization (fT2l ) in conjunction with the symmetry consequence ([TOl l 
leads to 

k{a,^)h{x,u;a) = k{-a,^)h{u,x;-a), 



hence 

h{x,u;cC) = h{u,x;—a). 
If we write the biquadratic h{x, u; a) as 

h{x,u,OL) = box^u^ + b[X^u + b2XU^ + b^x^ + b4U^ + bsxu + b^x + bju + b^, 

where bi = bi{a), then the coefficients obey the relations: 



(20) 



(21) 



hi 


-a) 


= ^2 (a) 


hi 


-a) 


= b4{a) 


b6( 


-a) 


= by (a) 


bo( 


-a) 


= bo{a) 


b5( 


-a) 


= bsia) 


b^[ 


-a) 


= b^ia). 


Equations ([15]) and (fT6l) give 






g{x,xuai 




g{x,xi;-ai,a2) 


g{x,xuai 


,a3) 


g{x,xi;-ai,a3) 



which implies 

g{x,xi;ai,a2) = cg{x,xi;-ai,a2) 
The factorization (fT2l) then leads to 

h{x,u;a) = ch{x,u;—a). 



(22) 
(23) 
(24) 



Comparing coefficients shows c = 1 and h is symmetric; that is, b2 = b\, b^ = b^, and 
bi = be. □ 



5 Derivation: Classification extension 

If the biquadratic h, given by ( fT9l ). depends on a it does so as an even function. However, 
if the only dependency on a is as an even function, we may consider a new dependency 
on a' := a^. In terms of a' (which is an 'undirected edge' lattice parameter) the new 
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symmetry proposed in this paper becomes the reflection symmetry of ABS, and leads to 
the ABS classification. Therefore, if the new symmetry is to lead to additional integrable 
lattice equations, augmenting the existing classification, then h must be independent of a. 

The classification of ABS ranks lattice equations according to canonical forms of the 
polynomial r{x) associated with them. To continue the parallels with the ABS classifi- 
cation (so that the place in the classification where the new lattice equations belong may 
easily be seen) we do the same. 

The simultaneous Mobius transformations x i— > {ax+b) / {cx + d) transforms the poly- 
nomial r as: 



Canonical forms of the quartic polynomial r{x) : 
r = 

r = I {r has one quadruple zero); 

r = x (r has one simple zero and one triple zero); 

r = x^ (r has two double zeros); 

r = x^ — I (r has two simple zeros and one double zero); 

r = 4x^ — g2X — g3, A = — 27^3 7^ (r has four simple zeros). 

From a given form of h, the associated lattice coefficients are found by solving the 
nine equations that arise by comparing coefficients of the biquadratic ([T2l ) (the lattice 
coefficients entering from the definition of g, dSjl) and the nine equations obtained by 
interchanging a and |3 and using (fT3] ) and the transformations of the lattice coefficients 
given in section|3] Considering the canonical forms for r(x) separately gives the following 
results: 



The only forms of h (that are independent of a) that give r = and lead to an inte- 
grable lattice equation obeying the assumptions of this paper are Mobius equivalent to 



The only forms of h (that are independent of a) that give r = x^ and lead to an in- 
tegrable lattice equation obeying the assumptions of this paper are Mobius equivalent to 
h = xu. 

Remaining canonical forms of r(x) 

Although there are biquadratics h independent of a that reduce to the particular r{x), 
there are no associated lattice equations obeying the symmetries given in section [3] 




(25) 



r = 



h = l. 
r = 
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5.1 Lattice equation of new symmetry associated with r = 

Solving 

g{x,u;a,^) = k{a,^)h{x,u;a) =k{a,^)h{x,u), 

g{x,u;^,a) =k{^,a)h{x,u;^) = -k{a,^)h{x,u), 

with h{x, u) = I, for the lattice coefficients, and imposing the symmetries given in section 
|3] leads to 

Q = K{a, p) (x - 3^) (v - u) +M(a, p) (« - v) + aM(-a, p) (x - 3;) +P{a, p) , (26) 

/t = -P(a,p)is:(a,p)-aM(a,p)M(-a,p), (27) 
and the conditions on the lattice coefficients: 



K{^,a) = -eK{a,^) , K{-a,^) = oK{a,^), 

M(p,a) = -8M(a,p) , M(-p,a) =£M(-a,p), 
P(P,a)=sP(a,p) , P(-a,p)=aP(a,p). 

This lattice equation is now a candidate for three-dimensional consistency. 

We now slightly re-write the lattice equation, and investigate the three-dimensional 
consistency condition, which will give us the final form of the integrable lattice equation. 
Dividing through by K{a,^) and defining m(a,P) := M{a,^)/K{a,^) and p{a,^) := 
P(a,p)//:(a,p), gives 

Q={x-y){v-u)+ m{a, p) (m - v) + m(-a, p) (x - y) + p{a, p) , (28) 

and the symmetry conditions become 



m(P,a) = wi(a,P) , m(-P,a) = -m(-a,P), 
p(p,a) = -7?(a,p) , p(-a,p) =7?(a,p). 

Investigating three-dimensional consistency gives two conditions that are necessary and 
sufficient for it to hold. The first is 

m(a3,ai) +ni(-ai,a2) -ni(a2,a3) = 0. (29) 

Using the symmetry conditions and differentiating ( [29l ) with respect to ttj gives 

9m(a3,ai) 3m(a3,a2) , 
da, - das 
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by a separation of variables argument. Further use of the symmetry conditions yields 

m{ai,aj) =f{ai)+f{aj). (30) 

The symmetry requirement m(— P, a) = — m(— a, P) shows /(— (x) = —/(a), hence, using 
the admissible transformations, we may take f{a) = a. Therefore m(a, P) = a + p. 

The second condition for three-dimensional consistency, after setting m{a, P) = a+ P, 
reads 

P(a3,ai)+p(ai,a2)+p(a2,cx3) =0. (31) 
Differentiating (|3TI ) with respect to a2 and using the symmetry conditions gives 

/7(ai,a2) =X(ai)-?i(a2). (32) 

The symmetries show X{a) may be any even function of a. (We may also replace X(a,) 
with a'^ where a- is a second lattice parameter associated with edge /.) Therefore, 

Q = {x-y){v - + (a + P)(m - v) + (p - a){x-y)+X{a) - ?i(p). (33) 

We call this lattice equation Bl. (Scaling such that the linear terms, (a + P)(m — v) and 
(P — a){x — y), disappear leaves the lattice parameters only appearing as even functions, 
hence the lattice equation becomes a member of the ABS classification, namely HI.) 

5.2 Lattice equation of new symmetry associated with r = x^ 

Solving 

g{x,u;a,^) =k{a,^)h{x,u;a) =k{a,^)h{x,u), 

g{x,u;^,a) = ^(P,a)/i(x,M;p) = —k{a,^)h{x,u), 

with h{x, u) = xu, for the lattice coefficients, and imposing the symmetries given in section 
[3] leads to 

Q = ea^:(-P, a)xu + e^:(P, a)vy + aK{-<x, ^)xv + K{a, ^)uy, (34) 

k = a^(-p, a)K{^, a) - aK{-a, p)^(a, p) , (35) 
and the condition on the lattice coefficients: 

K{^,-a)=aK{^,a). (36) 

This lattice equation is now a candidate for three-dimensional consistency. Investigating 
the three-dimensional consistency condition gives the two conditions: 

K{a2,ai) _ K{a2,a3) ^^^^ 



K{-a2,ai) K{-a2,(X3) 
(and its cyclic permutations) and 

K{ai,a3) K{a3,a2) ^(a2,ai) _ 
K{a3,ai) K{a2,a3) K{aua2) 



-8. (38) 



6 DISCUSSION 

A separation of variables argument from (|37] ) gives 

^(«2,ai) ^ fjai) 
K{-a2,ai) ~ f{-a2) 

Equation (1381 ) implies £ = — 1 and 

j^(ai,a2) _ A'(«i) 
/(:(a2,ai) ~ ^(a2)' 

So the symmetry condition (l36l) may be re-written as 

^ /(«2) ^ A'(«2) 

/(-a2) /j(-a2)' 

Hence 

2 = -A'(-P)jcM-A'(P)v3'+A'(-a)-^v+^(a)M3;. (42) 

If /J is an even function the lattice equation is part of the ABS classification (H3g=o)- If A* 
is not an even function we have an integrable lattice equation that is not part of the ABS 
classification. We call this lattice equation B2. 

6 Discussion 

That Bl and B2 are not part of the ABS classification is apparent by noting they are 
not of the ABS forms ^ and (01). Bl is a generalization of the lattice potential KdV. 
Setting ;u(a) = a in B2 it is the lattice potential mKdv. Setting ii{a) = —a — r in B2 it 
is the 'mixed lattice mKdV, discrete-time Toda lattice'. See lH for a survey of results on 
the lattice potential KdV, lattice potential mKdv, and 'mixed lattice mKdV, discrete-time 
Toda lattice' up to 1995. As already discussed, Bl can be reduced to HI and the form of 
B2 with an even function is H3g=o- Further historical remarks on the lattice equations 
in the ABS classification can be found in HI. 

Using transformations outside of our admissible, symmetry preserving, transforma- 
tions, B 1 can be transformed to HI and B2 to HSg^o; hence they should not be considered 
as new integrable lattice equations. (In the same way that, for instance, Al can be trans- 
formed to Ql.) Rather, the main result is that these are the only two integrable lattice 
equations with the given symmetry assumptions and directed edges. 

7 Conclusion 

The ABS classification, lU, has been extended by allowing the natural notion of directed 
edges. The full classification is shown in table [T] 

Our admissible transformations include simultaneous point change of all parameters 
that preserve the parity of the coefficients of the lattice equations as functions of the lattice 
parameters. This preserves the symmetry of the lattice equation. However, any simultane- 
ous point change of all parameters preserves the integrability, but destroys the symmetry. 
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(39) 



(40) 



(41) 
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HI 


Bl 


r(x) = 1 


Ql8=l 


Al8=l 


H2 




r(x) =x 


Q2 








r(x) = 


Q35=0 


A2 


H3 


B2 


r{x) =x^-8^ 


Q3 








r{x) = Ax^ - g2X - g3 


Q4 









Table 1: The classification of integrable lattice equations with symmetries 1, 2' and 3. 



This indicates that symmetry stipulations lose their importance on a deeper level. Table [T] 
is currently the largest list of integrable lattice equations; however, it would be desirable 
to have the full classification (or even the full list) of integrable lattice equations without 
symmetry assumptions. This remains an open problem. 
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